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Abstract 

We discuss invariant measures of partial differential equations such as the 2D Euler or Vlasov 
equations. For the 2D Euler equations, starting from the Liouville theorem, valid for iV-dimensional 
approximations of the dynamics, we define the microcanonical measure as a limit measure where N goes 
to infinity. When only the energy and enstrophy invariants are taken into account, we give an explicit 
computation to prove the following result: the microcanonical measure is actually a Young measure 
corresponding to the maximization of a mean-field entropy. We explain why this result remains true for 
more general microcanonical measures, when all the dynamical invariants are taken into account. We 
give an explicit proof that these microcanonical measures are invariant measures for the dynamics of 
the 2D Euler equations. We describe a more general set of invariant measures, and discuss briefly their 
stability and their consequence for the ergodicity of the 2D Euler equations. The extension of these 
results to the Vlasov equations is also discussed, together with a proof of the uniqueness of statistical 
equilibria, for Vlasov equations with repulsive convex potentials. 

Even if we consider, in this paper, invariant measures only for Hamiltonian equations, with no 
fluxes of conserved quantities, we think this work is an important step towards the description of 
non-equilibrium invariant measures with fluxes. 



*Author to whom any correspondence should be addressed. 



1 Introduction 



In a complex and chaotic dynamical system, such as a Hamiltonian system with a large number of degrees 
of freedom, or a non-equilibrium steady state where dissipation balances forcing on average, the knowledge 
of a non-trivial invariant measure is equivalent to the knowledge of the statistical properties of all physical 
quantities. It is thus an essential concept. In a turbulent problem, the knowledge of an invariant measure 
gives access to the stationary probability distribution function of all physical quantities, and gives a 
solution to the usual hierarchy of the n-point correlator of the velocity field, among other things. A series 
of very interesting recent works have proved the existence of invariant measures and described some of their 
properties, for instance in stochastic systems forced by noises [441 113 ISHl ES]- However, only in very few 
instances of complex systems, is an invariant measure explicitly known. Finite-dimensional Hamiltonian 
systems are among these exceptions: thanks to the Liouville theorem, a uniform measure on a constant- 
energy shell of phase space is invariant (microcanonical measure); the canonical Gibbs measures are other 
explicit examples of invariant measures. This essential remark is at the base of equilibrium statistical 
mechanics. 

For Hamiltonian partial differential equations, the situation is more complex. Indeed, the dimension of 
the system is then infinite. Thus, the meaning of phase space volume and the proper normalization of 
a uniform measure over a constant-energy shell of phase space are not clear notions. Microcanonical or 
Gibbs measures then have to be built carefully and their properties have to be checked. There are a few 
examples, where a Gibbs-type invariant measure has been proved to exist (see for instance |14| I15| in 
the case of the nonlinear Schrodinger equation, see also |45| I51j). As in |14| I15| . such a proof usually 
involves the study of approximations of the partial differential equations with finite dimension A^, and of 
limits of ensembles of measures when A^ goes to infinity. That the limit is actually an invariant measure 
of the initial partial differential equation completes the proof. We consider in this paper the construc- 
tion of microcanonical and Gibbs measures for the 2D Euler equations and for the Vlasov equations. We 
also consider other sets of invariant measures for these equations, and investigate their dynamical stability. 

The flow of a perfect fluid is described by the Euler equations, one of the oldest equations in mathemat- 
ical physics [30] . More than two and a half centuries after their discovery by Euler, these equations still 
offer great challenges to both mathematicians and physicists |30) . Two-dimensional flows and the two- 
dimensional Euler equations are mathematically much simpler than their three-dimensional counterparts, 
but still present some very interesting unsolved problems. The strong analogies between the 2D Euler 
equations and the Vlasov equations have been observed at least since the '50s: they are both nonlinear 
transport equations, the non-linearity being due to non-local and non-integrable (at long distance) inter- 
actions (long-range interactions) 0. We hope that this work may be very useful in the large number of 
systems with long range interactions that have been studied recently [9} \26 [ \TE [ I10 | [20]. 

One of the main physical phenomena arising from these two equations is the self-organization into large- 
scale coherent structures: large-scale particle clusters for the Vlasov equations with attractive potential, 
large-scale particle clouds (whose density proflles depend on both the dynamical invariants and the external 
confining potential) for the Vlasov equations with repulsive potential, or monopoles, dipoles, and parallel 

^Whereas in the 2D Euler and Vlasov equations, the long-range interactions and the associated non-additivity are essential 
for thermodynamical properties of the system (leading possibly to statistical ensemble inequivalence), they are not essential 
in the following discussion dealing with invariant measures and dynamics. For instance all the discussion of this paper can 
be easily adapted to the Quasi-Geostrophic model for which the interaction decays exponentially for distance much larger 
than the Rossby deformation radius. 
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flows for tlie 2D Euler equations. Such large-scale structures are analogous to geophysical cyclones, 
anticyclones, and jets in the oceans and atmospheres [13]. This analogy, understood thanks to the strong 
theoretical similarities between the 2D Euler equations on one hand, and the quasi-geostrophic or shallow- 
water models on the other hand, is one of the main motivations for the study of the 2D Euler equations. 
The 2D Euler equations also describe experimental flows, such as the transverse dynamics of electron 
plasma columns [65j, the large-Reynolds-number approximation of the dynamics of fluids when three- 
dimensional motion is constrained by a strong transverse field (rotation, transverse magnetic field in a 
liquid metal, etc. [B7J), or the dynamics of fluids in very thin geometries |58| . 

As first guessed by Onsager [57] . such a self-organization is naturally explained by equilibrium statistical 
mechanics: the infinite number of degrees of freedom involved should have a macroscopic behavior that 
corresponds to an average over a microcanonical measure. Because of the long-range interactions between 
vortices, such an equilibrium is not uniform and depends strongly on the boundary conditions. This 
explains the formation of vortices and jets. In his work [57], Onsager circumvented the difficulty of infinite 
dimensions discussed above, by studying the point- vortex model, which is a finite-dimensional Hamiltonian 
system, describing the dynamics of singular vortices. This model is actually a special class of solutions to 
the 2D Euler equations (please see [191 122[ |27] and references therein for recent developments of the kinetic 
theory of point- vortices) . The equilibrium statistical mechanics of the point-vortex model has a long and 
very interesting history, with wonderful pieces of mathematical achievements |57| [38] \T7\ l40l [28| \T9\ [3T| [32] . 

The study of the equilibrium statistical mechanics of the 2D Euler equations or of the Vlasov equations 
with usual smooth initial conditions implies considering the real infinite-dimension nature of the problem. 
An equilibrium statistical mechanics was proposed for these equations, based on the maximization of 
a Boltzmann-Gibbs entropy constrained by the conservation of the invariants, by Lynden-Bell for the 
Vlasov case [48] , and by Robert, Sommeria, and Miller for the 2D Euler case |59l [531 (see also |24) 
for a discussion on the analogies between both systems). Such an approach is basically a mean-field one: 
the use of the Boltzmann-Gibbs entropy for the entropy (or mean-field entropy) can be justified with 
the classical Boltzmann counting argument, which is meaningful only if correlations between the vorticity 
values (or values of the one-particle distribution in the Vlasov case) at different points can be neglected. 

Although this mean- field approach was a phenomenological assumption in the first papers |48| [59] [60] , 
it can be proved to be exact in systems with long-range interactions. For instance, physicists discussed the 
validity of the mean-field approach for the point-vortex model in the '80s, arguing that a Kramer-Moyal 
expansion when N goes to infinity is self-consistent. During the '90s, mathematicians and mathematical 
physicists proved that the mean-field approach is correct for the point-vortex model, using different tools 
\17\ [3T| I40| . For the 2D Euler and Vlasov equations, following the work of Bourgain for the nonlinear 
Schrodinger equation, such a proof should involve three steps: 

1. Proving Liouville theorems for finite-dimensional approximations of the dynamical equations, and 
defining microcanonical measures for those, 

2. Studying the infinite- limit for this set of measures and proving that the limit is actually described 
by the maximization of a mean field entropy, 

3. Proving that the limit measure is a dynamically invariant measure of the 2D Euler (or Vlasov) 
equations. 

For Galerkin truncations of the 2D Euler equations, e.g., using Fourier mode decomposition. Point 1 
above is a classical result [46] ; actually, the Euler equations verify a detailed Liouville theorem f52j. More 
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recently, Robert [6T] proved that a much larger class of approximations of the 2D Euler equations, obtained 
by L2 projections (e.g., finite-element approximations), verify a Liouville theorem. This last point is very 
important, because measures based on spatial truncations are much more natural than those based on 
Fourier mode decompositions for systems with long-range interactions. This makes the proof of Point 
2 much easier. Indeed, Michel and Robert |52) proved large-deviation results for ensembles of Young 
measures where the mean-field entropy appears as the opposite of the large-deviation rate function. A 
similar large-deviation result, justifying also the maximization of a mean-field entropy, was obtained from 
sets of measures based on spatial discretization [6l[7], but without any reference to the dynamics. These two 
results are essential steps, but give only a partial answer to Point 2, as the relation to dynamics is missing, 
making any further step towards Point 3 more difficult. As noted in |61| , another difficulty is that Fourier 
or L2 projections conserve only the quadratic invariants. Then, finite-dimensional measures that are based 
on the Liouville theorem and that take into account further invariants, which is the relevant procedure, 
are not invariant for the finite-dimensional approximation of the dynamical system. This makes Point 3 
out of reach. For instance, the derivation of Point 3 for the nonlinear Schrodinger equation by Bourgain 
[141115) relies deeply on the fact that this equation has only quadratic and linear invariants, and that finite- 
dimensional approximate measures are invariants for the finite-dimensional dynamical approximations. 

From this discussion, we conclude that although the validity of the mean-field approximation and the 
validity of the maximization of a mean-field entropy for some sets of natural measures have been proved 
[531 [SI E], a clear proof of the relation of those microcanonical measures with dynamics is still missing. 
Another important question is to know whether all measures built from any reasonable finite-dimensional 
truncation would lead to a unique limit measure in the limit of infinite dimension. 

As regards the relation between microcanonical (or canonical) measures and dynamics, we argue in this 
paper that it is not necessary to follow the classical program described above (points 1, 2, 3). More 
precisely, we argue that there is no logical need for the finite-dimensional approximate measures to be 
invariant measures of the finite-dimensional approximations of the dynamical system. We can indeed rely 
on the large-deviation results |53[ [6] and verify only a posteriori that the limit measure is an invariant 
measure of the 2D Euler equations. This is what we do in this paper, by studying directly the evolution 
equation for measures of the 2D Euler equations and studying ensembles of invariant measures. The 
discussion relies on Young measures, i.e., product measures for which the vorticity values at different 
points in space are independent random variables. The importance of Young measures in the context of 
equilibrium statistical mechanics was stressed in |52|,l61j. One of the main motivations for the present work 
is to give a dynamical meaning to Young measures, that we guess will be essential for further developments. 

We write the evolution equation for the characteristic functional of the vorticity field, for the 2D Euler dy- 
namics. Applying these equations to Young measures, we can study the class of Young measures that are 
invariant. We then note that microcanonical (or canonical) measures, which are actually Young measures, 
are dynamically invariant. The class of invariant Young measures is however larger; we explain that this 
property is related to a breaking of ergodicity for the 2D Euler equations. The dynamical approach of this 
paper also allows to discuss the dynamical stability of invariant measures, and particularly of invariant 
Young measures. The possibility to discuss dynamics is one of the main motivations for this work, as will 
be further discussed in the conclusion. 

In this paper, we also give a partial answer to the second question above: would any truncation lead to 
the same limit measure? For this we consider the case of the energy-enstrophy measure (which takes into 
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account only the energy and enstrophy as dynamical invariants). With direct explicit computations of the 
expectations for the Fourier A^-dimensional approximate measures, we prove that the limit measure is the 
same as the one obtained by maximizing the mean-field entropy (obtained as a large-deviation functional 
of a Young measure or from spatially discretized systems) . 

Energy-enstrophy measures were first investigated by Kraichnan |41) . in the canonical ensemble, and with- 
out considering the limit of an infinite number of degrees of freedom. In this paper, for the first time, we 
derive explicitly the microcanonical measure, which is the relevant physical one, for it takes into account 
the constraints explicitly. This derivation is based on expressing the constraint as a complex integral of 
an auxiliary variable, and analyzing precisely the integral asymptotics through the saddle-point approx- 
imation. This allows to actually take the limit of an infinite number of degrees of freedom with a fixed 
value of the invariants. This is a good opportunity to discuss again the Kraichnan theory: although the 
microcanonical measure leads to the same qualitative prediction as Kraichnan's, i.e., the condensation of 
energy in the largest-scale mode, we prove that microcanonical and canonical ensembles are not equivalent. 
More precisely, this is a situation of partial ensemble equivalence (see [29j for a definition). 

In section [21 we introduce the 2D Euler equations and their invariants, we define the microcanonical 
measure, and we give a pedagogic heuristic account of the meaning of the mean-field approximation and 
of the resulting statistical equilibrium distributions. 

In section [3l we perform the computations for the energy-enstrophy measure, both directly from its 
definition through Fourier modes and from the mean-field variational problem. We show that the results 
are equivalent and discuss ensemble inequivalence as well as the Kraichnan theory. 

Section m is devoted to the dynamics of measures for the 2D Euler equations. We characterize the set of 
invariant Young measures and their stability. In section [H we briefiy explain the generalization of these 
results to the Vlasov equations. We also prove the uniqueness of statistical equilibria in the case of convex 
repulsive interaction potentials. 

The relation of this work to non-equilibrium problems is discussed in section [U) 



2 Equilibrium statistical mechanics of two-dimensional flows 
2.1 2D Euler equations 

Defining the vorticity as w = (V x v) • , the 2D Euler equations take the simple form of a conservation 
law for the vorticity. They read 

(1) ^ + .-V uj = ; v = e^xV-0; uj = Aip, 

at 

where the solenoidal (incompressible) velocity v is expressed as the orthogonal of the streamfunction 
gradient Vip. We complement the equation uj = Aip with impenetrability boundary conditions: = on 
dD, where P is a simply connected domain. 

The kinetic energy of the fiow is conserved. It reads 

(2) S[^] = l [ dr v2 = i / dr (V^)' = [ dra;^. 
Other conservation laws are discussed in section 12.2.31 
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2.2 Microcanonical measures 



2.2.1 Theoretical foundations of equilibrium statistical mechanics 

Let us consider a canonical Hamiltonian system: {qi}i<i<N denote the generahzed coordinates, {pi}i<i<N 
their conjugate momenta, and H{{qi,pi}) the Hamiltonian. The variables {qi,Pi}i<i<N belong to a 2N- 
dimensional space called the phase space. Each point {{qi,Pi}) is called a microstate. The equilibrium 
statistical mechanics of such a canonical Hamiltonian system is based on the Liouville theorem, which 
states that the non-normalized measure 

i 

is dynamically invariant. The invariance of is equivalent to 
(3) E(I* + *)=». 



dqt dpi 



which is a direct consequence of the Hamilton equations of motion: 

dH 

OPi 

dH 

oqi 

Note that the equations of motion can also be written in a Poisson bracket form: 



(4) 



qi = {qi,H} , 
Pi = {pi,H} . 



Each term of the sum ([3]) actually vanishes independently: 

w. dqi dpi 
oqi opi 

Such a relation is called a detailed Liouville theorem. 

For all conserved quantities {Ii{q,p), ■ ■ ■ p)} of the Hamiltonian dynamics, the measures 

(5) iJc = ^Y{dpidqiF{Ii,...,Iri) 

t 

are also invariant measures, F being any function; Zc is a normalization constant. An important question 
is to know which of these measures are relevant for describing the statistics of the physical system. 

In the case of an isolated system, the dynamics is Hamiltonian and there is no exchange of energy or other 
conserved quantities with the environment. It is therefore natural to consider a measure that takes into 
account all these dynamical invariants as constraints. This justifies the definition of the microcanonical 
measure (for a given set of invariants {Ii{q,p), . . . , Iniq,p)})- 

1 " 

(6) fim {ho, ■■■,Ino) = 7777 j-^ n n ^ i^k{q,p) - 4o), 

Si {IIq, . . . , Ino) ■ j^^^ 
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where r2(/ig, . . . , J^q) is a normahzation constant — for smah {AI^}, . . . , InQ)YYk=i 

volume of the part of phase space such that for all k, Ik^ < Ik ^ Iko + ^Ik El- Then, the Boltzmann 

entropy of the Hamiltonian system is 

S = ks log , 

where fc^ is the well-known Boltzmann constant. When the system under consideration is not isolated, but 
coupled with an external bath of conserved quantities, other measures are necessary to describe properly 
the system by equilibrium statistical mechanics. A classical statistical mechanics result proves that the 
relevant functions in ([5]) are then exponentials (Boltzmann factors), and measures of type ([5]) are referred 
to as canonical measures (hence the 'c' subscript). The notion of a coupling with a single thermal bath (or 
to a bath of other conserved quantity than energy) in statistical mechanics assumes that if the system is 
coupled to a single bath, the total system (bath + system) is described by a microcanonical measure. This is 
an essential assumption in order to derive the canonical measure. By contrast, when systems are coupled 
to the environment through irreversible couplings (without detailed balance), or through several baths 
with different thermodynamic parameters, there is no reason anymore to expect a canonical measure to 
describe the statistics of the system. Fluid mechanics systems, where usually the stirring processes and the 
dissipation mechanisms are physical phenomena of a different nature, can never be considered as coupled 
to a single bath of some conserved quantities (in all works so far on fluid systems, where the distribution 
of energy or other conserved quantity have been studied, the distribution of conserved quantities is non 
Gaussian) . 

Hence, the relevant statistical ensemble for these models is the microcanonical one. In the following, 
we will work only with microcanonical measures as a base for all derivations. If equilibrium statistical 
mechanics is relevant for slightly non-equilibrium situations in fluid mechanics problems, it will most likely 
not be through canonical distributions, but through microcanonical distributions. 

In statistical mechanics studies, it is sometimes argued that, in the limit of an inflnite number of de- 
grees of freedom, canonical and microcanonical measures are equivalent. Thus, as canonical measures 
are more easily handled, they are preferred in many works. However, while the equivalence of canonical 
and microcanonical ensembles is very natural and usually true in systems with short-range interactions 
(commonly found in condensed matter theory), it is often not actually so in systems with long-range in- 
teractions, such as the 2D Euler equations (see for instance [9l 126 ^ \TE [ \TU [ [20 } l8l 129] and references therein). 

In statistical mechanics, a macrostate M is a set of microstates verifying some conditions. The conditions 
are usually chosen such that they describe conveniently the macroscopic behavior of the physical systems 
through a reduced number of variables. For instance, in a magnetic system, a macrostate M could be the 
ensemble of microstates with a given value of the total magnetization; in the case of a gas, a macrostate 
could be the ensemble of microstates corresponding to a given local density / (x, p) in the six dimensional 
space (x, p) (p, space), where / is deflned for instance through some coarse-graining. In our fluid problem, 
an interesting macrostate will be the local probability distribution p (x, a) da to observe vorticity values 
CO (x) = (J at x with precision do". 

more natural definition of tlie microcanonical measure would be the uniform measure on the submanifold defined by 
Ik = Ikg for all k. This would require adding determinants in formula and imply further technical difficulties. In most 
cases though, in the limit of large number A'^ of degrees of freedom, these two definitions of the microcanonical measure 
become equivalent. Indeed, the measures have then large-deviation properties (saddle-point evaluations), where N is the 
large parameter, and such determinants become irrelevant. We note that in the original works of Boltzmann and Gibbs, the 
microcanonical measure refers to a measure where only the energy constraint is considered. 
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If we identify the macrostate M with the values of the constraints that define it, we can define the probabil- 
ity of a macrostate P (M) dM. If the microstates are distributed according to the microcanonical measure, 
P (M) is proportional to the volume of the subset of phase space where microstates {(?i,Pi}i<j<^ re- 
alize the state M. The Boltzmann entropy of a macrostate M is then defined to be proportional to the 
logarithm of the phase space volume of the subset of all microstates {^i,Pj}x<j<Ar that realize the 
state M. 

In systems with a large number of degrees of freedom, it is customary to observe that the probability of 
some macrostates is concentrated close to a unique macrostate. There exist also cases where the probability 
of macrostates concentrates close to larger set of macrostates (see for instance [39j). Such a concentration 
is a very important information about the macroscopic behavior of the system. The aim of statistical 
physics is then to identify the physically relevant macrostates, and to determine their probability and 
where this probability is concentrated. This is the program we will follow in the next sections, for the 2D 
Euler equations. 

In the preceding discussion, we have explained that the microcanonical measure is a natural invariant 
measure with given values of the invariants. An important issue is to know if this measure describes also 
the statistics of the temporal averages of the Hamiltonian system. This issue, called ergodicity will be 
discussed in section 4.4. 

2.2.2 Hamiltonian structure for the 2D Euler equations 

The first step to define the microcanonical measure is to identify the equivalent of a Liouville theorem 
and of the dynamical invariants. The Euler equations describe a conservative dynamics. They can be 
derived from a least-action principle |64| I36| . just like canonical Hamiltonian systems. It is thus natural 
to expect a Hamiltonian structure for them as well. There are however fundamental differences between 
infinite-dimensional systems like the 2D Euler equations and canonical Hamiltonian systems: 

1. The Euler equations consist in a dynamical system of infinite dimension. The notion of volume of 
an infinite-dimensional space is meaningless. Hence, the microcanonical measure cannot be defined 
straightforwardly. 

2. For such infinite-dimensional systems, we cannot in general find a canonical structure (pair of canon- 
ically conjugated variables {(??,, Pi} describing all degrees of freedom). It exists however a Poisson 
structure: one can define a Poisson bracket {•,•}, like in canonical Hamiltonian systems (jl]), so that 
the dynamics writes 



where % is the Hamiltonian. 

For infinite-dimensional Hamiltonian systems (such as the 2D Euler equations), the Poisson bracket in 
([7]) is often degenerate |35^ 155] . leading to the existence of an infinite number of conserved quantities. 
These conservation laws have very important dynamical consequences, as explained in the next section. 
A detailed description of the Hamiltonian structure of infinite-dimensional systems is beyond the scope of 
this paper. We refer to |351 155| for a description of the Poisson and Hamiltonian structure for many fluid 
systems. In the next two sections, the dynamical invariants and the Liouville theorem are discussed in the 
context of the 2D Euler equations. 



(7) 
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2.2.3 Casimir conservation laws 



2D Euler equations ([T]) conserve an infinite number of functionals, named Casimir invariants (or Casimirs 
for short). They are all functionals of the form 

(8) CM= I drs(L^), 

Jv 

where s is any sufficiently smooth function. As mentioned in section 12.2.11 Casimir conserved quantities 
are related to the degeneracy of the Poisson structure in infinite-dimensional Hamiltonian systems. They 
can also be understood as the invariants arising from Noether's theorem, as a consequence of the relabeling 
symmetry of fluid mechanics (see for instance [64]). 

Let us note A{a) the area of T) with vorticity values less than a, and 7 [a) the vorticity distribution: 

(9) ^ " ^'^^^^ Iv ' 

where XB is the characteristic function of the set B C "D, and |D| is the area of D. Since equations 
([1]) express transport by an incompressible flow, the area 7 (a) occupied by a given vorticity level a (or 
equivalently A{a)) is a dynamical invariant. 

The conservation of the distribution 7 (a) is equivalent to the conservation of all Casimir functionals ([8]) . 
The domain-averaged vorticity Q, enstrophy Q2, and higher moments of the vorticity {Gk}k>3 are Casimirs 
of particular interest: 

(10) for A; > 1, Gk [u] = [ drw^ (with G [uj] := Gi [u]). 

Jv 

Note that if T) is bounded, G is also the circulation: G = Jg-p v • dl . 

In any Hamiltonian system, symmetries are associated with conservation laws, as a consequence of 
Noether's theorem. Then, if the domain T) is invariant under rotations or translations, there will be 
conservation of angular momentum or linear momentum respectively. If the domain displays such sym- 
metries, these conservation laws have to be taken into account in a statistical mechanics analysis. 



2.2.4 Detailed Liouville theorem and microcanonical measure for the dynamics of conser- 
vative flows 

In order to discuss the detailed Liouville theorem, and build the microcanonical measure, we decompose 
the vorticity field on the eigenmodes of the Laplacian on T>. We could decompose the field in any other 
orthonormal basis. The Laplacian and Fourier bases prove simpler for the following discussion, whereas 
finite-element bases are much more natural to justify a mean-field approximation and to obtain large- 
deviation results for the measures, as will be discussed in section 12.41 

We call {ej}i>i the orthonormal family of eigenfunctions of the Laplacian on T): 

(11) - Aej = XiCi , / dr CjCj = 5ij . 

Jv 
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The eigenvalues {Aj} are arranged in increasing order. For instance, for a doubly periodic or infinite 
domain, {e,} are just the Fourier modes. Any function g defined on T) can be decomposed into g{t,r) = 
Ei5i(*)ej(r) with gi{t) = J^dr g{t,r)ei{r). Then, 

+00 

a;(r,t) = ^uji{t)ei{r). 

i=l 

From ([H), 

(12) UJi = AijkiOjUJk , 

where the explicit expression of Aijk need not be known for the following discussion. For (I12p . a detailed 
Liouville theorem holds: 

(13) V^, ^ = 

(see |46| . [42]). Note that even though we have discussed the detailed Liouville theorem in the context of 
mode decomposition, more general results exist [61)|72| H. 



Microcanonical measure From the detailed Liouville theorem (jl3p . we can define the microcanonical 
measure. First, let us define the n- moment microcanonical measure: 

1 " 

(14) {E,T,,..., r„) = — n ^ [^\ [Qk M - Ffc) , 

where 8 is the energy ([2]) and \Tk} are the vorticity moments (jlOp . the subscript ra still standing for 
'microcanonical'. A precise definition of //m,n requires the definition of approximate finite-dimensional 
measures: for any observable (j)M depending on M components {ijJi}i<i<M of w, we define 

^ „ Af n 

{fJ'Z,n^ 4>m) = TTT-p fTT W^UiS {£n [w] - -E) JJ 5 {Gk^N [uj] - Ffc) (/)M (^^1, • • • , ^m) , 

^ ^n,N [-t^ ■, 1 J ■ ■ ■ 1 n ) J ^_-|^ 

where and {Gk,N} are finite-dimensional approximations of lE" ([2]) and {Gk} (|10P respectively. Then we 
define {lJ,m,n,4'M) = ^^^N^ooifJ-m m 4'm) ■ As explained in the next paragraph, Qn,N has usually no finite 
limit when N goes to infinity, so the definition of the normalization factor il„ (E, Fi, . . . , r„) in the formal 
notation (jl4p implies a proper rescaling. 

{Mm,n}n>i IS expected to be a set of invariant measures for the 2D Euler equations (this is easily verified 
through formal computations). The microcanonical measure corresponding to the infinite set of invariants 
{^k}k>i then defined by 

(E, {Ffe}) = hm 

n— f oo 

and is denoted by 

^ oo oo 

(15) {E, {Tk}) = ^^^^^^^^^ n d'^i ^ t'^] - ^) n ^ ['^i - ta.) • 



direct consequence of the detailed Liouville theorem (|13p is that any truncation of the Euler equations also verifies a 
Liouville theorem |42| . This result is actually much more general: any approximation of the Euler equations obtained by an 
L2-projection on a finite-dimensional basis verifies a Liouville theorem (see [61]). For truncations preserving the Hamiltonian 
structure and a finite number of Casimir invariants, see 1721. 
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Equilibrium Boltzmann entropy The normalization factors [E, Ti, . . . , r„) and Q, {E, {Ffc}) define 
the Boltzmann entropies 

(16) Sn{E,ri, . . . ,Tn) = kBlognn{E,Ti, . . . ,Tn) and S {E,{Tk}) = kBlogn{E,{Tk}) . 

The behavior of expressions like J Hi^i 6 {£n [w] — E) 13^=1 ^ {Gk,N ['^] ~ Tfc), for large A^, is expected 
to be typically of the form C{N) exp [NS^ {E, Fi, . . . , F^)). In the definition of r2„, the prefactor C{N) is 
omitted (the entropy is defined up to a constant independent of the physical variables), so Sn and S are 
actually 'specific entropies' (entropies per degree of freedom). 

2.3 Validity of a mean-field approach to the microcanonical measures 

In the previous section, we defined the microcanonical measure for the 2D Euler equations. In this section, 
we give a heuristic explanation of the reason why a mean-field description of the microcanonical measure is 
exact and give references for more precise results (large deviations for sets of measures). At the core of our 
discussion about mean-field approaches, lies the result that if the microstates are distributed according to 
the microcanonical measure, the probability distributions of vorticity at different points are independent 
(product measure). This is not only a result of precise mathematical works (large deviations for sets of 
measures), but also the deep reason for the validity of the mean- field approach. This is also the main 
reason for the interest of Young measures for the 2D Euler equations |52| . 

Because the large deviations of sets of measures are rather technical results from probability theory, for 
pedagogic reasons, we study the energy-enstrophy microcanonical ensemble in section IXTl with elementary 
mathematical tools. We prove in Appendix D that the correlation coefficient between Ci;(r) and t<;(r') is 
zero for the energy-enstrophy microcanonical measure. It could be proved without much difficulty that, in 
addition, u;(r) and u){r') are actually independent variables. As said above, the statistical independence 
of vorticity values at different points is a much more general result and is essential. Let us first analyze an 
extremely important implication: the possibility to quantify the phase space volume (Boltzmann entropy) 
through the Boltzmann-Gibbs formula. 



Boltzmann entropy of a macrostate and Boltzmann Gibbs formula A classical example where 
degrees of freedom can be considered independent is an ensemble of particles (say, hard spheres) undergoing 
collisions in the dilute limit (Boltzmann-Grad limit |69|). Microscopically, particles travel with typical 
velocity v and collide with each other after traveling a typical distance I, called the mean free path. Let a 
be the diffusion cross-section for these collisions. One has a = ira^, where parameter o is of the order of the 
particle radius. The mean free path is defined as I = l/(7ro^n), where n is the typical particle density. The 
Boltzmann equation applies when the ratio a/Z is small (Boltzmann-Grad limit). In the limit a/ 1 — )• 0, 
any two colliding particles can be considered independent (uncorrelated) as they come from very distant 
regions. This is the base of Boltzmann's hypothesis of molecular chaos (Stosszahl Ansatz). It explains why 
the evolution of the ^u-space distribution function /(x, p, t) may be described by an autonomous equation, 
the Boltzmann equation (the ^u-space is the six-dimensional space of spatial variable x and momentum p) . 
In statistical mechanics, a macrostate M is a set of microstates verifying some conditions. The conditions 
are usually chosen such that they describe conveniently the macroscopic behavior of the physical system 
through a reduced number of variables. The Boltzmann entropy of a macrostate M is defined to be 
proportional to the logarithm of the phase space volume of the subset 0,m of all microstates {%,Pi}i<j<Ar 
that realize the state M. In the case of a dilute gas, the distribution / (x, p) can be identified with the 
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macrostate: it is the set of aU possible microstates {xi,Pi}^^-^j^ such that the number of particles in the 
volume element AxAp around (x, p) is / (x, p) AxAp (a precise mathematical definition goes through 
the limit N ^ oo, see [33]). We note that the Boltzmann entropy of the subset of phase space with fixed 
invariants is the equilibrium Boltzmann entropy defined in section 12.2.41 formula p6|) . 

There is a classical argument by Boltzmann (which can be found in any good textbook on statistical 
mechanics) to prove that the Boltzmann entropy of the distribution / is, up to a multiplicative constant, 
given by the Boltzmann-Gibbs formula: 

(17) S[f] = - I dxdp/log/. 

We stress that this formula for the Boltzmann entropy is not a Gibbs entropy @. The essential point 
is that formula (jl7p is a valid counting of the volume of the accessible part of phase space, only when 
particles can be considered independent. For instance, for particles with short-range interactions studied 
by Boltzmann, this is valid only in the Boltzmann-Grad limit. 

As discussed above, for the microcanonical measure of the 2D Euler equations, vorticity field values are 
independent. As we will explain below, the reason is completely different from the Boltzmann case: there 
is now no dilute-gas (Boltzmann-Grad) limit. Nevertheless, the consequences will be the same: if we 
define p (r, a) such that p (r, a) drdo" be the probability to have values of u! between a and a + dcr in the 
area element dr around r, then the entropy 

r r+co 

(18) S[p] = - dr da pin/), 

Jv J~oo 

actually quantifies the phase space volume. In order to give a precise meaning of this last sentence, we 
first define the mean-field microcanonical variational problem. 



Mean-field microcanonical variational problem As p is a local probability, it verifies a local nor- 
malization 

/ + 00 
dap{a,r) = l. 
-oo 

The average vorticity, for probability density p, is 

/ + 00 
da ap (cj, r) . 
-oo 

The average vorticity (|20p is related to the average streamfunction ip so that w = Atp. 

The conservation of all Casimir functionals jH]), or equivalently of the known vorticity distribution ([9]), 
imposes a constraint on the local probability density p: 

(21) D[p]{a)^ [ drp(a,r) = 7(a). 

Jv 

*The Gibbs entropy S = —k J p{pi, Qi) log2(p(pi, Qi)) dptdqt is an ensemble entropy, a weight on the phase space, whereas 
the Boltzmann-Gibbs entropy is an integral over the /i-space. In the case of dilute gases, the Boltzmann-Gibbs entropy is just 
the opposite of the H function of Boltzmann. We avoid this terminology here since our discussion is not related to rela:xation 
towards equilibrium, and because the equivalent of an H theorem has never been proved for the 2D Euler equations. 
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Then the mean-field entropy S of the system is given by the variational problem 
(MVP) S{E, 7) = sup {S[p] \£[0] = E,D[p] = j}. 

{p|7V[p]=l} 

where £ [a)] is the energy ([2]) of the average vorticity field ui. 

An essential point has to be noted about the energy constraint in (jMVPp : the constraint is expressed in 
terms of the average vorticity field (hence the expression, 'mean-field approximation'), meaning that cor- 
relations between vorticity values at different points are negligible, and meaning also that when computing 
the energy, fluctuations around the average may be neglected. 

An essential point is that, up to addition of constant terms (i.e., independent of the physical parameters), 
the mean-field entropy (jMVPp is exactly the same as the Boltzmann entropy defined from the rescaled log- 
arithm of the phase space volume, in equation (I16p . The definition of the entropy (fTBl) and the variational 
problem (jMVPp seem so different, that the fact that they express the same concept is astonishing. This 
type of results is indeed one of the great achievements of statistical mechanics. In section 13.11 we show 
that it is verified in the case of the energy-enstrophy measure, using explicit elementary computations. 

Why is the mean-field entropy equal to the Boltzmann entropy? The deep reason why vorticity 
field values are independent for microcanonical measures, and henceforth why entropy can be expressed 
by (134p can be explained rather easily at a heuristic level. Correlations between variables could appear 
through the dynamical constraints only: energy, Casimirs, and so on. For instance, the energy of the 2D 
Euler system can be expressed in a form where interactions between vorticity values appear explicitly, 
using the Laplacian Green function H (r, r') (defined by AH (r, r') = 6{r, r') with Dirichlet boundary 
conditions), we have: 

(22) £[uj] = -- [ druj{A-^uj) = -- [ [ drdr' u (r) H{r,r')uj{r'). 

In formula (j22p above, H (r, r') appears as the coupling between vorticity at point r and vorticity at point 
r'. The Laplacian Green function in a two-dimensional space is logarithmic, hence non-local. Thus, a;(r) 
is coupled to the vorticity at any other point of the domain, not only close points. 

For people trained in statistical mechanics, it is natural in systems where degrees of freedom are coupled 
to many others, to consider these degrees of freedom statistically independent at leading order, and a 
mean-field approach should be a valid approximation. For example, in systems with nearest-neighbor 
interactions, a mean-field approach becomes exact in high dimensions, when the effective number of 
degrees of freedom to which one degree of freedom is coupled becomes infinite. For people not trained 
in statistical mechanics, this can be understood simply: when the number of coupled degrees of freedom 
increases to infinity, the interaction felt by one degree of freedom is no more sensitive to the fluctuations 
of the others, but just to their average value, owing to an effect similar to what happens for the law of 
large numbers. Then a mean-field treatment becomes exact, which is equivalent to saying that different 
degrees of freedom may be considered statistically independent. 

Because of the non-locality of the Green function, the vorticity field at one point is virtually coupled to an 
infinite number of degrees of freedom, and then a mean-field treatment is exact. This also explains why 
the energy appearing in the variational problem (jMVPp is computed from the average vorticity field. 
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To formalize the preceding heuristic explanation, in order to prove that the mean-field approximation is 
exact and that the Boltzmann-Gibbs formula (jl7p is relevant, we need a rather technical discussion. We 
will not explain this in detail. This was justified by theoretical physicists for the point-vortex model in 
the seventies (assumed to be valid by Joyce and Montgomery [38J and later proved to be self-consistent in 
a Kramer-Moyal expansion). In the eighties, rigorous mathematical proofs were given also for the point- 
vortex model (see |3 H HO l [T7] and references therein). In the modern formulation of statistical mechanics, 
the entropy appears as a large-deviation rate function for an ensemble of measures, justifying (fT7|) and 
the variational problem (jMVPp . The proof of such large-deviation results leading to the microcanonical 
measure for the 2D Euler equations, justifying the mean-field approach, can be found in j52j (see also [6] 
and references therein). 



We thus conclude that a mean-field approach to the microcanonical measure of the 2D Euler equations 
is valid. This justifies the use of entropy (jlSp and of the variational problem (|MVPp . This step is a 
crucial one as it leads to a drastic simplification compared to a direct computation from the definition of 
the microcanonical measure (jl4p -(jl5 p . The first presentation of the equilibrium statistical mechanics of 
the 2D Euler equations in this form dates from the beginning of the '90s with the works of Robert and 
Sommeria, and those of Miller |59^ [5^ [60] 162] . Thus, we call this theory the Robert-Sommeria — Miller 
(RSM) theory. 



2.4 Solutions to the mean-field variational problem for the microcanonical measure 

The aim of this section is to describe the critical points of the mean-field variational problem (IMVPp . 
following the first papers \59\ [5^ l60l I62| . For this purpose, we use the Lagrange multiplier rule to take 
account of the constraints: the first variations of (IMVPp . 

r r+oo 

6S - drA{r)6N{r) - p6£ - / daa{a)6j{a) = 

are zero for any perturbation 5p, where A[r), /3, and a{a) are the Lagrange multipliers associated with 
the conservation of A^(r), £ and 7(cr) respectively. We obtain that the probability density distribution p 
verifies the Gibbs state equation: 

(23) p{a,r) = - with (u) = / da exp (irn - a (a)) . 

[p^ (r)) 7_oo 

We see that p depends on r through the average streamfunction '\\) only. From (I20p and ()23p . we see that 
there is a functional relation between the equilibrium average vorticity and the streamfunction: 

(24) w = 5 (/3V5) with g{u) = ^ log Zo,{u). 

This last equation characterizes the statistical equilibrium. It should be solved for any value of (/3,a((T)). 
Then, one has to compute the energy and vorticity distributions as functions of /3 and a(cj). For given 
energy E and distribution 7(0"), among all possible values of (/3, a((T), p((t)) solving (|23ll24p . the maximizer 
of the entropy (jMVPp is selected. 
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3 Energy-enstrophy microcanonical measure for the 2D Euler equations 



The energy-enstrophy microcanonical measure is defined as 

^ oo 

(25) ^J.m,K (E, Ts) = ^Y^^Y^ n '^(^ [^] - M - r2). 

This is the measure where only the quadratic invariants are taken into account. There is a priori no 
physical reason to exclude the other invariants; however, the energy-enstrophy microcanonical measure can 
be interesting, because it is, in some cases, a good approximation of the complete microcanonical measure. 
Our real motivation to treat it in detail is rather pedagogical: it will be very useful to prove with this simple 
example, using elementary explicit computation, the equivalence between the microcanonical measure 
introduced in section r2.2.4l through Fourier mode decomposition, and the solution to the microcanonical 
mean-field variational problem of section 12.41 

We compute the entropy and the probability distribution function for the amplitude of each mode. These 
computations are performed, always in the microcanonical ensemble, on one hand directly from the defini- 
tion of the energy-enstrophy measure (section 13. ip , and on the other hand from the mean- field variational 
problem (section l3.2p . 

The energy-enstrophy measure was treated and discussed at length by many authors in the '70s, including 
Kraichnan (see [42j, a precise discussion can be found in [49]). However, these computations were always 
performed in the canonical ensemble. The following discussion gives the first derivation in the microcanon- 
ical ensemble, and the first observation of ensemble inequivalence for the energy-enstrophy ensembles — 
microcanonical and canonical (section 13. 3p . The energy-enstrophy ensembles are an elementary example 
of the so-called partial equivalence [29], in the theory of ensemble inequivalence. We shall come back to 
discuss Kraichnan-type results in section [331 



3.1 Direct computation of the energy-enstrophy measure from its finite-dimensional 
approximation 

Following the discussion of section 12. 2. 4^ the energy-enstrophy microcanonical measure is defined through 
A^-dimensional approximations: 

1 ^ 

(26) Hm,K = lim n^ j^ with /i^^ = TT duji 6{Sn [^] - E)d{g2,N [w] - Fa), 

N-^oo IIk^N [E,L2) fj^ 

where we use the same notation as in section [2.2.4[ so 28^ [w] = "^iLi^f / K ^^nd F2,Ar [oo] = "^j^i^of. In 
the following, we assume that the first mode is non-degenerate: Ai 7^ A2 (this is always true for simply 
connected bounded Lipschitz domains, but this is wrong for doubly periodic boundary conditions in a 
square domain). 

The main technical difficulty is to compute 

N 

(27) nK,N {E, Fa) = / n ^i^N [co] - E)6{g2,N [co] - F2), 

1=1 
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and the entropy 
(28) 



Sk (E, T2) = lim - log [Qk,n {E, T2)] - C{N, {AJ), 

Af— s-oo iV 



where C does not depend on the physical parameters. It depends only on and on the geometric factors 
{Ai}, and can be discarded as the entropy is always defined up to an arbitrary constant. 

The computation of ^k,n and Sk, using representation of the delta function as an integral in the complex 
plane, is given in Appendix B. It yields the result 



VLk,n{E,T2) ~ C^{N,{\i])C^{{\i}j2,N) 

N~¥oo 



exp [NSk{E,T2] 



^iog(r2-2Aiii;) + i^ 



(29) with Sk {E, r2 

(see ()73p and (|74p . page l36p . where C3 (A^, {Aj}) does not depend on the energy or enstrophy, and C4 has 
no exponentially large contribution (lim jv -5.00 (logC4) /N = 0). 

We now describe finite- A'^ effects for //^^ (I26p . the finite- A'^ approximation of the energy-enstrophy mi- 
crocanonical measure ^m,K (|25p . It is easy to see from (j26p . that for the A^-dimensional measure /U^j^, 
the distribution function for the amplitude Un of mode e„ is given by 

^K,N-i-M {E - ul/2X„.,r2 - ul) 



(30) 



^k,n{E,T2) 

where the definition of i^K,N-i;X„ is the same as that of ^k,n (|27p . but with integration over Un excluded, 
and with constraint < max {2An-E, r2}. 

The distribution function for the energy En = Li;^/2A„ of mode is obtained through the change of 
variable -P/v.„ {En) dEn = PN,n i^n) dcj„. Using result (j29p for both i^K,N-i;Xi (then Ai has to be replaced 
with A2) and ^k,n, we obtain 

(31) P.AEr) ~ crexp[Ariog(r2- 2A2i. + 2(A2-A,)i.,)/2] ^ < < 

and Pn,i (Ei) = otherwise. C is a normalization constant which does not depend on Ei. From this 
expression, we can see that the most probable energy is Ei = E. Moreover, the distribution is exponentially 
peaked around Ei = E, so that in the infinite- A/^ limit (energy-enstrophy microcanonical distribution), we 
have 

Pi{E,) = 6{E-Ei). 

This is a striking result: for the energy-enstrophy microcanonical measure, all the energy condensates in 
the first mode. 



The type of result (j3ip is called a large-deviation result. It describes accurately the distribution function 
for El, however large the deviations from the most probable value, in the limit A^ — t- 00. If large deviations 
in E — El are disregarded, a good approximation for large A^ of the finite- A^ distribution is the exponential 
distribution: 



(32) (^1 



N^oo 



exp 



-A^ 



r2~2XiE 

\/ E — El 



{E-Ei 



for Q<Ei<E and N^/^ {Ei -E)^l. 
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We note that the distribution for coi is exponential as weh. The amphtude of the departure of Ei from 
value E is thus proportional to 1/N and to — 2XiE) / (A2 — Ai) . 

The distribution of the energy £"„ of mode e„ is obtained similarly as 

(33) P.AEn) - c^exp [ATlog (T, - 2X,E- 2 (A„ - E.)] ^ ^ ^ ^ 

N-^oc VEn 

For infinite A^, the energy-enstrophy microcanonical distribution is thus a delta function with zero energy: 

PniEn) = S{En). 

Disregarding large deviations, the finite-A^ distribution is also well approximated by an exponential dis- 
tribution (this time a Gaussian distribution for a;„), with typical departure from of order 1/N for the 
energy and a variance of order 1/y/N for ujn- It may also be checked that for large n (A„ ^> Ai), the 
variance of the enstrophy becomes independent of n (asymptotic equipartition of the enstrophy). 

Results such as (f3T]) - (j33p are classical in statistical mechanics: typical departures from the most probable 
value have a Gaussian distribution with variance of order 1 /^/N, except for variables whose most probable 
values are at the edge of the accessible range. In the latter case (for instance, En or oji in the example 
discussed above), the distribution is exponential with typical departure of order 

We note that for the 2D Euler equations, only the infinite- limit is relevant, and finite- effects have 
no dynamical counterpart. They may be of interest for truncated systems only. 



From the preceding discussion, we see that all the energy is concentrated in the first mode, and that the 
excess enstrophy r2 — 2Ai-E' goes to smaller an smaller scales, leading to zero energy and zero enstrophy 
in every mode except the first one. This condensation of energy in the first mode is the main physical 
prediction of the microcanonical energy-enstrophy ensemble. 



3.2 Energy-enstrophy microcanonical measure from a mean-field approach 

Let us compute the entropy in the (microcanonical) energy-enstrophy ensemble, now starting from the 
mean-field variational problem, and compare the results with those of section 13.11 

The mean-field variational problem in the (microcanonical) energy-enstrophy ensemble is the equivalent 
of (jMVPp but with only quadratic invariants taken into account: 

(34) Sk {E, T2) = sup I -^S[p] \£[oJ] = E , [ drda p = j 

{p\N[p] = l} U^l J J 

Note that we seek here to maximize the specific entropy because this is what leads to the 

actual measure of phase space volume. It is customary in the literature to ignore the 'i/\D\ prefactor for 
convenience, as done for instance in (jMVPp . 

In order to compute the critical points of the constrained variational problem (j34p . we introduce Lagrange 
multipliers j4(r),/3, and a, associated with the conservation of N{r)^£., and Q2 respectively (see (I19p . ([2]) 
and (jlOp . respectively, for the expression of these quantities). Critical points of (j34p are such that 



5S_ 

W\ 



[ dr A{r)6N{r) - ^58 - aSQi = V 5/9. 
Jv 
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This is equivalent to 

p{a,r) = p*{a,r) = S(r)el^l(^"'^W-""'), 
where the prefactor B{r) = exp(— 1 — j4(r)|D|) is determined from the normaUzation constraint: 



(35) = y^e-l-l{'-^)'. 

The computation above, yielding a Gaussian distribution in the energy-enstrophy ensemble, is a classical 
result noted in many previous works (see, for instance, [23j). 

Substituting expression (j35p into (jlSp . we get 

(36) ^5[p*] = -^loga 

for the expression of the entropy (see calculation in Appendix B-2). We thus conclude that the maximum- 
entropy solution will be the one verifying the constraints with minimum value for a. 

We now compute E and r2 as functions of /3 and a. For this, we compute the average vorticity using ([35|) 
in 1^: 



r+co g _ 

uj{t)= dcr(Tp*((T,r) = — VXr) 



We thus have a;(r) = /\ip{r) = /3'i/'(r)/(2a), so we can deduce that vorticity and streamfunction are 
proportional to a Laplacian eigenmode: for some n > 1, u){r) = Acni^) ™d ^(r) = — ^/Anen(r), with 
I3n = -2a„A„. From £[u)\ = E, we have E = ^V(2A„) and ip{r) = -^lEjXn e„(r). From r2 = 
Jj^^ d(Ta"^p*(a", r) (see Appendix B-2 for detailed computation), we get r2 = 2A„i? + l/(2a„). Thus, we 
can see that the minimum value a„* of a is obtained for n* = 1. 

Finally, since entropy (|36p is maximum for a > minimum, the first eigenmode (n = 1) is the one selected. 

We are left with ^ 

tD(r) = \l 2\\E ei(r) and a 



2(r2-2Ai^) ' 
so that the equilibrium entropy is 

Sk {E, T2) = \ log (r2 - 2\^E) + . 

Comparing this result with (j29p . we can conclude that the entropy computed from the mean-field vari- 
ational problem, in the energy-enstrophy ensemble, is the same as the one computed directly from the 
definition of the energy-enstrophy microcanonical measure, through Fourier mode decomposition (finding 
a)(r) = y/2y[Eei{r) is equivalent to finding Pi{Ei) = 6{E - Ei) and Pn{En) = 5{En)). 



3.3 Ensemble inequivalence 

From the entropy, we can compute the inverse temperature /3 = dSx/dE = — Ai/ (r2 — 2Ai£') < 
and fugacity a = dSK/dT2 = 1/ [2 (r2 — 2Aii?)]. These thermodynamical coefficients are related through 
/3 = — 2Aia. This relation shows that some couples of thermodynamical coefficients are not obtained in the 
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energy-enstrophy microcanonical ensemble, in contrast to what would be expected in the thermodynamics 
of classical condensed matter systems. Moreover, the determinant of the Hessian of Sk, that is, d'^Sx/dE'^- 
d'^Sx/dr'^ — {d'^ Sk /dEdT2)'^ , is zero, showing that Sk is not strictly concave, unlike what would be 
expected for an entropy in the case of short-range interacting systems. Both these properties are signs 
of non-equivalence between the microcanonical and canonical ensembles: the two ensembles would give 
different predictions (see for instance |9l |29l [26]). This case of ensemble inequivalence, for the energy- 
enstrophy ensembles, is actually a case of partial equivalence (see |29) for a definition). 

A detailed discussion of ensemble inequivalence and related phase transitions, for statistical equilibrium 
with linear relation between vorticity and streamfunctions, including the case of the energy-enstrophy 
ensemble just discussed can be found in |70) . 

3.4 Comments on the Kraichnan energy-enstrophy theory 

The term 'condensation' was proposed by Kraichnan from the analysis of the energy-enstrophy canonical 
ensembles |42j . As explained in section [2.2.11 canonical measures are not relevant for fluid systems; they 
may be useful only when yielding results equivalent to the ones from microcanonical measures. Kraichnan 
noticed this and worked nonetheless with canonical ensembles, maybe because he did not know how 
to perform microcanonical computations, most likely because at that time the possibility of ensemble 
inequivalence was nearly unknown H. Unfortunately, as mentioned earlier, the energy-enstrophy ensembles 
display an instance of partial ensemble inequivalence. These remarks explain the difficulties encountered 
by Kraichnan when analyzing the canonical measures, and why he wrongly concluded that a statistical 
mechanics approach would work only for truncated systems. Working in microcanonical ensembles actually 
allows to build invariant measures for the real (non-truncated) 2D Euler equations. If one were interested 
in truncated systems, then Kraichnan's work would remain very useful. 

More importantly, when looking closely at Kraichnan's works (see, for instance, |42) page 565), one sees 
that in the canonical ensemble, a complete condensation of the energy on the gravest mode occurs only for 
specific values of the thermodynamical parameters. For most values of the thermodynamical parameters, 
an important part of the energy remains on the other modes. Still Kraichnan argued, probably from 
numerical observations available at the time and from physical insight, that these cases leading to a 
condensation were the most interesting ones. The microcanonical treatment we propose here proves that a 
complete condensation occurs whatever the values of the energy and of the enstrophy, in the microcanonical 
ensemble. A complete condensation is actually observed in many numerical simulations. We thus conclude 
that the physical insight of Kraichnan and his concept of condensation describes the relevant physical 
mechanism, but that a treatment in the microcanonical ensemble provides a much better understanding, 
and overcomes the preceding contradictions. 

Limitations of the energy-enstrophy approach There is a priori no reason to consider only the 
energy and enstrophy invariants, except for being able to solve the mathematics easily. From the discussion 
of section (j2.4p . we know that a mean-field approach is exact for the microcanonical measures in the case 
of 2D Euler equations, and thus the description of any microcanonical measure (corresponding to any set 
of invariants) is not difficult. 

^The first observation of ensemble inequivalence was made in the astrophysical context [481 134) . while a thorough study 
[261 [9] and understanding of the importance of ensemble inequivalence for two-dimensional flows [661 1401 1311 1291 170| are more 
recent. 
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As has been shown in previous works (see, for instance, |12)). when taking into account ah invariants, the 
energy wih no longer be hmited to the first mode ei. The energy-enstrophy measure may still be a good 
approximation in some cases: in the limit of small energy, for instance, most of the energy will remain 
in the first few modes. The notion of condensation will thus remain valid only roughly speaking, at a 
qualitative level. 

By contrast, in some cases such as that of doubly periodic domains with aspect ratio close to but different 
than one (see |12j). the notion of condensation would lead to completely wrong predictions. 

4 Invariant measures of the 2D Euler equations 

In the previous section, we have built microcanonical measures for the 2D Euler equations and argued 
that they are a special set of Young measures. In this section, we consider the dynamics of measures and 
more specifically the dynamics of Young measures. We give a direct proof that sets of Young measures, 
including microcanonical ones, are invariant measures of the 2D Euler equations. For this we derive the 
evolution equation verified by the characteristic functional, in section 14.11 This equation describes the 
evolution of all the statistics of the system: it is equivalent to the Liouville equation and includes the 
hierarchy of equations describing the evolution of the statistics of the vorticity field. 

This section is completely independent from the previous one. The proof that sets of Young measures are 
invariant is independent from the building of the microcanonical measure in the previous section. However, 
the construction of microcanonical measures from the Liouville theorem, in the previous section, shows 
that among the set of Young measures, microcanonical measures have a specific meaning. 

4.1 Evolution of the characteristic and cumulant-generating functionals 



cumulant-generating function h{l) = log/(/), where the angle brackets denote average over the measure 
of the random variable x. In order to describe the temporal evolution of the statistics of the vorticity 
field uj{r,t), it will prove very useful to use a generalization of the characteristic and cumulant-generating 
functions to random fields. We consider an ensemble of initial conditions {ujq (r)}. Each of these initial 
conditions evolves according to the 2D Euler equations, defining an ensemble of solutions of the 2D Euler 
equations {uj{r,t)}. We define the characteristic and cumulant-generating functionals of the ensemble 
{loq (r)}, respectively, as 



where the angle brackets denote ensemble average over realizations of cjq. We define similarly the charac- 
teristic and cumulant-generating functionals of the ensemble {uj{r,t)}, respectively, as 



where the angle brackets still denote ensemble average over realizations of the initial conditions ujq. 

We now use that each realization uj{r,t) is a solution to the 2D Euler equations ([T]), in order to derive the 
evolution equation for F. A straightforward computation, reproduced in Appendix D-1, leads to 



For any random variable x, it is customary to define 






(37) 
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where G is the Green function for the velocity: 

(38) v(r) = J dr' G(r, r')a;(r'). 

We note that the evolution equation for the characteristic functional ([37|) is a linear equation, as is the 
classical Liouville equation. 

The equation for the cumulant-generating functional is also obtained straightforwardly (see Appendix 
D-2): 

p'" -u+'JJ ^'<'' ■ °<'' ' ' [amm ^W)W) 

4.2 Young measures and their dynamics 
4.2.1 Young measures 

We recall that Young measures are uncountable product measures: the probability distribution of the 
vorticity field at an arbitrary number of points {r^} is given by the product of the independent measures 
p{a,Yk) at each point r^. We note that at each point, p is normalized 



/ + 00 
do- p(cj,r) = 1. 
-oo 



As we see below, the fact that vorticity values at different points are independent variables has important 
consequences. 

The set of deterministic vorticity fields are oj (r) is a special class of Young measures with p (a, r) = 
5{a-u:{v)). 

The set of microcanonical measures described in section (j2.4p . is a special class of Young measures. They 
are defined as 



(40) p^,{c,}(o-,r) 
where Zo,{u) = J^^dfi e'^"-"^'^). 



1 



Z{p^ (r)) 



If the local probabilities pi((T, r) and p2{(T,r) define two Young measures, then for any < / < 1, Ipi + 
{1 — l)p2 defines a Young measure. 



4.2.2 Cumulant-generating functionals for Young measures 

Let us evaluate the cumulant-generating functional of a Young measure. For this purpose, we consider 
h (/,r), the cumulant-generating function of the local probability p{cr,r) at each point r: 

/ + 00 
da e^'Xa,r), 
-oo 
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and the average vorticity field 

/'+°° dh 
(42) ^(r) = j ^ daap{a,r) = —{0,T). 

Let us now consider a sufficiently regular A(r) (for instance, Riemann-integrable for the purpose of the 
following discussion). The characteristic functional of the entire vorticity field is easily computed using an 
approximation by a finite Riemann sum: 

F[X] = (e'^-f ^^"^'^^"^'^"'^ = lim (e''^^k=iHrk)'-^{rk)'^ ^ 

where |P| still denotes the domain area. Since variables {^(r/j)}^ are independent, we use the fact that 
the characteristic function of a set of independent variables is equal to the product of the single-variable 
characteristic functions: 

N 



k=l 

SO that 



I — 1 ^ / 



H[X] = log F[X] = ^hm f; h (^A(r,), r,) = lun^ ^ E | (C r.) A(r,) , 

k=i ^ ^ k=i 

where we have used /i(0, r^) = 0. Then, using (I42p . we have 

(43) H[X] = i J A(r)w(r)dr. 

The cumulant-generating functional is linear in A. Actually, comparing (I43p to (j47p . we see that for 
sufficiently regular A(r), the cumulant-generating functional is the same as that of a deterministic field 
with vorticity the average vorticity ci;(r). Hence, the statistics of any observable that is obtained as the 
domain integral of a sufficiently smooth function of a;(r) depends only on the average vorticity 0{r). This 
can be seen as an example of a law of large numbers for an infinite sum of independent variables. 

In the following, we call an observable S smooth additive, if S is obtained as an integral over the vorticity 
field: 

(44) S = J dr(^(r)w(r). 

Then, from the preceding result, we can conclude that the distribution of any smooth additive observable 
S is a delta function. As an illustration, we compute the characteristic functional of the velocity field at 
point r, using A(r') = AG (r, r'), where G is the velocity Green function (j38p . Then, 

(45) log ^e'^^^"")^ = H[XG (r, r')] = iX J d'r G (r, r') 0{r') = iAv (r) . 

We see that the cumulant-generating function of v (r) is linear in A. The velocity field has no fiuctuations, 
and hence is a delta function centered at the average value v (r). Once again, this is nothing but the law 
of large numbers. 
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We can easily generalize the computation of the cumulant-generating functional (j43p to classes of non- 
regular A(r). As an example, we consider A (r) = L (r) + 15 {r — tq), where L (r) is Riemann-integrable 
and I is a scalar. A direct generalization of the computations preceding (|43p leads to 

H[L{r) + W(r - ro)] = i J Lir)u{r)dr + h (/, tq) . 

Using this result with L{r) = X(f), we can describe the joint probability of any smooth additive variable S 
(|44p and of the vorticity at point tq. We conclude that S and u (rg) are independent random variables, 
u) (ro) having distribution p (a, r) and S having a delta distribution centered at the average value S = 
J dr (p{r)Lj{r). 

This result can be extended to the joint probability distribution of the vorticity at an arbitrary number of 
points and of any smooth additive observable. For instance, we can conclude that for Young measures, the 
velocity field is a random variable independent of the vorticity field, having a delta distribution centered 
at the average velocity v(r), with Ai/j = O. 



4.2.3 Dynamics of Young measures 



We obtain the dynamics of /i(A,r), the cumulant-generating function of the vorticity field at point r 
(|4ip . either from (j39p or by direct averaging of the 2D Euler equations. Using that the velocity field is 
independent of the vorticity field, as explained at the end of section |4. 2. 2( we obtain 

dh 

(46) + V • V/i = . 

ot 

From this equation, it is clear that any initial measure which is Young measure remains a Young measure 
over time. Moreover equation (j46p describe the whole dynamics, as for Young measures, any observable 
can be derived from h. The dynamics of Young measures is thus rather simple. 

For the dynamics of Young measures, we could have worked directly with p. The evolution equation for p 
is just 

dp „ 

However it is more convenient to work with H (and h) as soon as perturbations to a Young measure are 
considered. 



4.3 Classes of invariant measures 
4.3.1 Deterministic dynamical equilibria 

For any dynamical system, equilibria of the deterministic equations are trivial invariant measures. This is 
obviously also the case for the 2D Euler equations. 

Let us consider a stationary solution of the 2D Euler equations ti;o(r), and the associated velocity field 
vo(r), such that vo(r) ■'Vujo{r) = 0. The associated characteristic and cumulant-generating functionals are 
Fo[A] = eiI<^^Mr)^o{r) 

(47) Ho[X] = i [ dr \{r)ujo{r) 
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respectively. It is easily verified that Fq and Hq are equilibria of (j37p and (|39p . respectively, as expected. 
Indeed, the second term in the l.h.s. of equation (|39p is then 

-i J drVA(r) • vo(r)a;o(r), 

which is trivially null. This is also trivially checked from (j46p . 

4.3.2 Invariant Young measures 



Prom (|39p and the cumulant-generating functional (j43p . we see that a necessary condition for a Young 
measure to be invariant for the 2D Euler equations is that for any sufficiently regular A(r), 



J drVA(r) • v(r)w(r) =0. 



Thus, a necessary condition is that the average vorticity ti)(r) be a dynamical equilibrium of the 2D Euler 
equations. From this equation, we see that h is transported by the average velocity v. Then a further 
necessary condition for a Young measure to be invariant is that h be invariant over any streamline of the 
velocity v. 

This is also a sufficient condition. Indeed, in 2D Euler equations, vorticity is just transported by the 
velocity field. Then, because for a Young measure, the velocity has no fluctuations, if the velocity is 
moreover stationary and if the distribution does not depend on the streamline, then the Young measure 
is invariant. Then, invariant Young measures are the ones for which h is invariant over streamlines of v. 

A smaller class of invariant Young measures of interest is the one for which p depends in a functional way 
on the streamfunction ip = ip : p = p{a,ip (r)). This property has to be self-consistent: 



CO = Alp = J da ap{a,ilj (r)). 



4.3.3 Microcanonical measures 



The set of microcanonical measures (I40p is a special class of Young measures. From the previous compu- 
tations, their cumulant-generating functional is 

A(r)Li;(r)dr with io{r) = g (r)) and g (u) = — log Z . 

Because of the functional relation between vorticity and streamfunction to'(r) = g (^(r)), w(r) is a dy- 
namical equilibrium of the 2D Euler equations. In conclusion, the microcanonical measures are invariant 
measures of the 2D Euler equations. 

4.3.4 Quasi-invariant Young measures 

We now consider the class of Young measures for which Ci)(r) is a dynamical invariant of the 2D Euler 
equations, but for which /i(A, r) is not invariant over each streamline. Then, from (j46p . because velocity 
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v(r) is stationary, /i(A, r) is just transported along each streamline. Therefore, from a microscopic point 
of view, such a Young measure is not invariant, but from a macroscopic point of view, it is: any smooth 
additive observable of the vorticity field, including the velocity field, is invariant. We call such a measure 
a quasi-invariant Young measure. 

4.4 Ergodicity 

Section [2] describes the statistical equilibria through the variational problem (jMVPp . The solution of 
this variational problem is the most probable state and also, thanks to the large-deviation property, the 
state around which an overwhelming majority of states do concentrate, for the microcanonical measure. 
Besides, the microcanonical measure is the most natural invariant measure of the 2D Euler equations with 
the dynamical constraints. 

Having described a natural invariant measure of the equations is an important theoretical step. Another 
important point would be to know if this invariant measure is the only one having the right values 
for the dynamical invariants. The evolution of one trajectory of the dynamical system also defines a 
measure (through time averaging). If we knew the invariant measure were unique, then it would mean 
that averaging over the microcanonical measure is equivalent to averaging over time. When this uniqueness 
property holds, we call the dynamical system ergodic. 

Generally speaking, the ergodicity of a dynamical system is a property that is usually extremely difficult 
to prove. Such proofs exist only for very few extremely simple systems. Ergodicity is actually thought to 
be wrong in general. For instance, in Hamiltonian systems with a finite number of degrees of freedom, 
there often exist islands in phase space in which trajectories are trapped. The common belief in the 
statistical mechanics community is that those parts of phase space where the motion is trapped exist, but 
occupy an extremely small relative volume of the phase space, for generic systems with a large number of 
degrees of freedom. Apart from a few systems which were proved to be integrable, this common wisdom 
has successfully passed empirical tests of a century of statistical mechanics studies. 

There is no reason to suspect that this general picture should be different in the case of the 2D Euler 
equations, in general. It is thus thought that an overwhelming number of initial conditions will have 
a dynamics consistent with the microcanonical measure predictions. However, similarly to most other 
Hamiltonian systems, the 2D Euler equations are actually non-ergodic, the proofs being extremely simple, 
given the discussion on invariant Young measures in the previous section. 

Indeed, in section 14.11 we have proved that any Young measure for which oj (r) is a stationary solution of 
the 2D Euler equations is either an invariant or a quasi-invariant measure. The class of invariant measures 
corresponding to ensemble of trajectories with given values of the invariants, is then much, much larger 
than the class of statistical equilibrium invariant measures with the same invariants. This proves that 
nontrivial sets of vorticity fields are dynamically invariant. In this restricted sense, this proves that the 
2D Euler equations are not ergodic. 

This theoretical argument proving non-ergodicity is in accordance with previous remarks about the phe- 
nomenology of the 2D Euler or quasi-geostrophic equations. For instance, it was observed numerically 
that initial conditions with localized vorticity, in large domains, remain localized (see |23] and references 
therein; [23] actually proposes an interesting phenomenological modification of the microcanonical measure 
approach to cope with this localized dynamics problem). Another example of possible non-ergodicity is 
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the dynamics close to stable dynamical equilibria of the equations. When trajectories come close to such 
equilibria, they can be trapped (frozen) as was seen in some numerical simulations. A classical argument 
by Isichenko |37| is that for initial conditions close to parallel flows, 'displacement in certain directions 
is uniformly small, implying that decaying Vlasov and 2D fluid turbulence are not ergodic'. Even if the 
predicted algebraic laws by Isichenko are most probably wrong, the fact that displacement in directions nor- 
mal to the streamlines is uniformly small is probably right, thus being another argument for non-ergodicity. 

An important point to be noted, is that the Navier-Stokes equation with stochastic forces can be proved 
to be ergodic |16| . This ergodicity refers to invariant measures of the Navier-Stokes equations, which 
are non-equilibrium invariant measures with fluxes of conserved quantity. A very important point is to 
understand the limit of weak forces and dissipation for such invariant measures and to study their relations 
with the invariant measures of the 2D Euler equations. Some very interesting results can be found in [44) . 

4.5 Stability of invariant measures 

As discussed in the previous section, the dynamical stability of invariant measures is an essential point. In 
the following two sections, we give a first discussion of the dynamical stability of invariant Young measures. 

We first discuss, in section [4.5.11 the stability of invariant Young measures when the perturbation is such 
that the initial condition is still a Young measure. We conclude that the stability then depends only on 
the unperturbed average vorticity ljq (r) and velocity vo(r). More precisely, if in the framework of the 2D 
Euler equations small perturbations of a) (r) lead to finite Lagrangian transport in the direction transverse 
to the streamlines, then the Young measure is stable to perturbations among Young measures. This 
condition of finite Lagrangian transport cross to the streamlines, is true for the 2D linearized equations 
for a whole class of parallel flows, including flows with stationary streamlines |11) . Moreover, we guess 
it is also true for the nonlinear 2D Euler equations, for large classes of parallel flows and stable vortices, 
even if there are still no proofs yet. 

In section I4.5.2| we consider the stability of dynamical equilibria of the 2D Euler equations, subjected 
to perturbations whose measure is not necessarily in the the class of Young measures. For this study, 
we assume perturbations to be small. We prove that the linear stability of the unperturbed flow in the 
framework of the 2D Euler equations, implies the linear stability in a statistical sense. 

4.5.1 Stability among Young measures 

As explained in section 14.2.31 the dynamics of a Young measure is equivalent to the dynamics of its 
cumulant-generating function h{X, r) (|46p . where v(r) is the velocity field corresponding to a)(r). Moreover, 
the set of Young measures is dynamically stable. We can thus consider stability of Young measures among 
the set of Young measures (we take the initial condition as a Young measure that is close to the unperturbed 
invariant one /io(A,r)). We note that this includes the case when each realization leading to the Young 
measure is perturbed by the same initial perturbation 6uj (r), in which case /i(A,r) = /io(A,r) + iX6uj (r). 

We have seen that invariant Young measures are the ones for which v(r) is stationary and h{X, r) is 
constant over every streamline of v(r). We consider such an invariant Young measure, denoting the 
associated cumulant-generating function by /io(A,r) and the associated velocity by vo(r), and a small 
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perturbation to : h = + 5h. We recall (see (j46p ). that the dynamics is 



dh 

(48) _+v.V/i = 0, 

with u){v) = 1^ (0, r) and v(r) the corresponding velocity field. From ()48p . we have 

(49) ^ + V • = 0. 

at 

We conclude that u) satisfies the deterministic 2D Euler equations. A necessary condition for /iq to be 
stable is thus that wq be stable for the 2D Euler equations. 

A variety of notions of stability exist for the 2D Euler equations, depending on the norms used to control 
the initial conditions and the evolving solutions. As far as the initial conditions are concerned, from a 
physical point of view, it is very natural to consider small perturbations of the initial velocity, as small- 
scale vorticity fiuctuations are usually not controlled. We note that (I48p and (l49l) are readily solved using 
Lagrangian coordinates: defining R(t) by 

dR 

= v(R) with R(t = 0)r = r. 

at 



Then, using the incompressibility of v, we deduce thatR(t) is invertible at each time t, and we have: 
w(r, t) = u}{n{t)-^Y, 0) and h{\, r, t) = h{\, R(t)"^r, 0). 

As discussed in paragraph 14.3.41 the motion along streamlines does not matter for the invariance of the 
measure. A natural definition of stability is therefore imposing that the motion normal to the streamlines 
remain small over time. We call this Lagrangian stability; it can be defined more precisely as, say. 



Ve, 3a: ]- j dr {5wf {t = < a || V^^wq • R(t)r||< e, 

where V^wq is the unit vector orthogonal to Va;o and where different norms ||.|| define different notions 
of stability. 

Then, clearly, from the preceding discussion, because the dynamics of Young measures is just the transport 
by the average velocity, and h{X, r, t) = h{X, I{.{t)~^r, 0), it is natural to define the Lagrangian stability of 
the Young measure, just as to be equivalent to the Lagrangian stability of the of 2D Euler equations for 
the average vorticity field. 



As can be readily seen from the results in ^llj . a whole class of parallel flows, including flows with 
non-monotonic velocity profiles, are Lagrangian-stable, as far as the linear dynamics is concerned (if 
the perturbed velocity field evolves according to the linearized 2D Euler equations, then the associated 
Lagrangian transport in the direction transverse to the streamlines is uniformly bounded over time and 
proportional to the initial perturbation amplitude). Even if there is no proof yet, we guess that this is 
also true for the (nonlinear) 2D Euler dynamics for a large class of parallel fiows or stable vortices. This 
would prove the stability of Young measures for the nonlinear dynamics. 

Then, any further discussion needs a detailed study of the relaxation (asymptotic stability) of the 2D 
Euler equations, which is not available yet but will be considered in future works. 
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4.5.2 Stability of invariant measures to small velocity perturbations 



We consider now the stability of a deterministic solution to the 2D Euler equations and the effect of small 
perturbations. The statistics is not limited to Young measures. 

We start with the 2D Euler equation which is verified for any realization: 



dt 



+ V • VJ] = 0, 



where il. = + euj , V = Vq + ev, and Vq • Vf^o = (e is the perturbation amplitude). Developing these 
equations yields 



(50) 



duj 
'dt 



+ L[uj] + ev • Vo; = 0, 



where the operator L, giving the linearized 2D Euler equations, is defined by 



(51) 



L[u] = Vo • Vw + V • Vno = Vo • + dr' G(r, r')w(r')^ • VJIq 



The evolution equation for the characteristic functional F[X\ = (^e^ I is then 



(52) ^+11 dr'dr VA(r) 



G{r,r'){ei 



6^F 



6X{r)6X{r') 



6F 



5F 



0. 



(a detailed derivation is provided in Appendix D). Likewise, the cumulant-generating functional H = lnF 
satisfies 



(53) 
dH 
'dt 



+ // dr'dr VA(r) 



G(r,r') 



ei 



5^H 



+ 



6H 6H 



5X{r)6X{r') 6X{r) 5X{r') 



6H 

6X{?': 



Vo(r) 



6H 
6X{r) 



0. 



We expand the equation for the cumulant-generating functional in powers of e: H = HQ + eHi+e^H2 + ■ ■ ■■ 
At lowest order (e°) we have 



(54) 

where the linear operator C is defined by 
" 6H 



dt 



+ C 



5Hn 



6X{r) 



0, 



6X{r) 



(55) 



y"drVA(r)- Vo(r)^ + J]o(r)y"dr' G(r,r' 
dr A(r)L 



JH_ 

6X{r') 



6H 
aw 



(see equation (jSip ). We now remark that Hq satisfies the same equation as the one for the cumulant- 
generating functional in the case of linearized 2D Euler equations. 
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At the linear level (first order in e), the stability of the measure is thus equivalent to the stability of 
the operator £. We remark that a detailed knowledge of the properties of L is sufficient to describe the 
properties of £; indeed, 



As a consequence, the stability of L (exp {tL) uniformly bounded over time) is a necessary and sufficient 
condition for the stability of C (exp uniformly bounded over time). In the case of stable parallel 
base flow Vq = U (y) e^, a detailed study of the asymptotic behavior of the linear operator L is provided 
in in particular, for any perturbation, it is proved that ||exp [tL) \\ is usually composed of the 

contribution of few modes (often no mode, actually) plus the contribution of a continuous spectrum that 
decays algebraically for large times. These results hold in the cases of both monotonic and non-monotonic 
velocity profiles [/, and are probably easily generalizable to the case of stable circular vortices. 

We note that the expansion of H at higher orders could be performed easily. Discussion of the convergence 
of such an expansion requires a detailed study of the relaxation (asymptotic stability) of the (nonlinear) 
2D Euler equations, which is not available yet and that will be considered in future works. 

As already stated several times, we note that a complete theory for the stability of Young measures 
requires the understanding of the nonlinear relaxation of the 2D Euler equations. Even at a linear level, 
a more complete study of the stability of Young measures would involve investigating the effect of small 
perturbations on any invariant Young measure, not only on deterministic solutions as done in this section. 
This requires a more involved expansion than the one in this section and will also be considered later on. 

5 Invariant measures of the Vlasov equation 

In this section, we consider the Vlasov equation. For the sake of simplicity we consider one-dimensional 
physical systems: however, the discussion easily extends to any dimension. Thanks to the the theoretical 
similarity between the 2D Euler and the Vlasov equations, noted decades ago, all the discussions about 
the microcanonical measures. Young measures, invariant measures, and stability of invariant measures of 
sections [21 13. H and 14.11 easily extend to the case of the Vlasov equation. In the following sections, we 
only describe briefly the statistical equilibrium measures and the invariant Young measures. We also give 
a proof of the uniqueness of statistical equilibria in the case of repulsive convex potentials, which can be 
useful in many future studies. 

5.1 The Vlasov equation 

We consider a set of particles subjected to their mutual two-body interactions with potential W . Each 
particle located at point x is subject to the potential (/'discrete (2;) = X^i^i W{x — Xj), where {xi} are the 
particle positions. W is an even function. Classical physical arguments and mathematical proofs justify 
that, when W is regular enough, it is natural to consider the following continuum approximation to this 
potential: 




(56) 
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The time evolution for the one-particle phase space distribution function f{x,p,t) satisfies the Vlasov 
equation, given by 

(57) ^+p^_^^ = 

dt dx dx dp 

If {xi} were N independent random variables distributed according to the distribution /, equation ([56|) 
would then follow from the law of large numbers and would be a good approximation to i?5>discrete up 
to corrections of order \/^/N. Replacing the true discrete potential by (j) thus amounts to neglecting 
correlations between particles (the equivalent of the Stosszahl Ansatz) and finite- effects. The potential 
'/'discrete being replaced by an average one, namely (j), may be seen as a mean-field approximation to the 
dynamics. 

As can be easily verified, the Vlasov equation (I57p inherits the conservation laws of the microscopic 
Hamiltonian dynamics, for instance, for the energy 



(58) H[f]= / dxdp 



^2 2 

and for the linear or angular momentum, when the system has the corresponding translational or rotational 
symmetry, respectively. 

If we define ip[x,p) = —{p'^/2 + </>(x)) and 

dil^ d'4>\ d V f d d 
dp^ dx J \dx' dp 

then, the Vlasov equation ([57|) can be recast into 

(59) ^ + v.V/ = 0, 

with V • V = (see equations ([1])). This simple remark explains the deep analogy between the 2D Euler 
and the Vlasov equations. Like the vorticity for 2D flows, / is transported by an incompressible flow. This 
explains most of the following properties. 

The Casimir functionals 

(60) Csif] = J dxdps{f{x,p,t)) 
are invariant for any function s. 

Let fm = sup{/}. We denote A{a) the area of D with / values greater than a, and 7(0") the vorticity 
distribution: 

dA - f 

(61) ^^^^"~da = /p^^^{'^</(^)</-} 

(see section [2.2.31 for the deflnition of xb) H- The area 7(0") of a given /-level a (or equivalently A{a)) is 
a dynamical invariant. The invariance of 7 (a) is equivalent to the invariance of all Casimirs Cs[/]. 



^These definitions are different from those of section [2.2.31 (2D Euler equations), because with the Vlasov equation, the 
area of D is infinite. 
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For the dynamical equilibria, there is no time dependence for (j) and ^. Then from (j59p we conclude that 
any distribution for which the distribution functions are constant over isovalues lines of "i/; are dynamical 
equilibria. For instance, for any /q, distribution of the type f{x,p) = /q {ip{x,p)) are dynamical equilibria 
(we note that this relation has to be self-consistent as the potential defining (j) is also computed from /). 

5.2 Equilibrium statistical mechanics of the Vlasov equation 

The equations for the equilibrium statistical mechanics of the Vlasov equation were first written by Lynden- 
Bell |48| for self-gravitating systems; extensions and discussions of the analogy with the 2D Euler equations 
for the equilibrium statistical mechanics were first discussed in |24| . The equilibrium statistical mechanics 
can be considered following exactly the same steps as for the 2D Euler equations, in section [2j Then a 
mean field approach will be valid and the equilibrium distribution will be a Young measure characterized 
by p{a,x,p) the probability distribution for / to take the value a at the point {x,p) of the phase space, 
with the normalization da p{a, x,p) = 1. 

Let us define the average one-particle distribution function 

/+00 
dcr ap{a,x,p). 
-oo 

The mean-field variational problem defining the statistical equilibrium is 

(62) SiE,^) = sup \s[p] I H[f] = E, D[p]ia) = [ dxdpp{a,x,p) = 7(a) | 

{Pl I^^ da p{a,x,p) = l} I J } 

(see (|MVPP ). where = — J dxdp plogp. The equilibrium probability density distribution then reads: 
where Z„(u) = da e"^-^^") . 

5.3 Invariant measures of the Vlasov equation 

Let us now look at the evolution of the characteristic functional, defined by 

p]^y^ _ l^^i S dxdp\{x,p)f{x,p)^ 

Using a Green function formalism, in order to keep the analogy with the 2D Euler equations, we have 

v[/](x,]5) = j dx'dp'Gix,x';p,p)f{x',p) + {0,p), 



where 



G(x,x';p,p')=(o,-^(x-x')). 
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By analogy with (j37p . we get 



— + J dxdp X{x,p)p- 



6F 



5X{x,p) 



+i I dxdpdx'dp — 
op 



dW 
dx 



(x — x) 



6^F 



6X{x,p)6X{x' ,p') 



+ 



0. 



The dynamics of Young measures, characterized either by p {a, x,p), or h (A, x,p) = log (J da p {a, x,p) exp (i 
is given by 

+ v. V/i = 0. 

ot 

The set of invariant Young measures is thus the set of Young measures such that h (or p) is constant over 
any isoline of the average particle energy ip. We can define quasi-invariant Young measures similarly to 
the case of the 2D Euler equations, as discussed in section 14. 3.41 



As in the case of the 2D Euler equations, because the set of Young measures and hence the set of invariant 
measures is much larger than the set of equilibrium measures, the Vlasov equation is non-ergodic. 

The discussion of the ergodicity in the framework of the Vlasov equation has a long history, starting with 
the works of Lynden-Bell |48| . A lot of recent works have made detailed comparisons of the prediction 
of the equilibrium statistical mechanics with numerical simulations [3l IH [21 [H [TH i see also a detailed 
discussion in [21j and references therein. The qualitative results are similar to the ones for the 2D Euler 
equations: whereas some cases definitely show not so good a prediction due to the equilibrium statistical 
mechanics, because of lack of ergodicity, in most cases equilibrium statistical mechanics provides a fairly 
good prediction of the final self-organized state. For instance, this theory has been used to predict the final 
bunching parameter of a free-electron laser [5]. Moreover, the prediction skill of the equilibrium theory is 
expected to be better when the dimension increases. We stress that the Vlasov equation in unbounded 
physical space without a confining potential presents some specific difficulties, as then the microcanonical 
measure is not defined. 



6 Perspectives 

In this paper, we have studied classes of Young measures which are invariant measures of the 2D Euler 
equations. These classes include microcanonical and canonical equilibrium measures, but not only. Our 
approach was to consider the problem directly from a dynamical perspective, by looking at the evolution 
equations for the characteristic functional and for the cumulant-generating functional. 

Our main motivation and the interest of this approach is, first, to study the stability of invariant measures 
and, second, to be able to generalize the results to other dynamical systems, for instance the 2D Navier- 
Stokes equation with stochastic forces. 

In sections 14.5.11 and I4.5.2| we began the study of the stability of invariant Young measures. We have also 
stressed that a more complete study of their stability requires new results about the relaxation towards 
dynamical equilibria of the 2D Euler equations. Very interesting results have been obtained recently for 
the relaxation towards dynamical equilibria of the Vlasov equation [56] . We hope that similar results will 
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be obtained soon for the 2D Euler equations, and that these wih be used in the future to complete the 
study of the stabihty of both the 2D Euler and Vlasov equations. 

The class of invariant measures we describe in this paper are no-fluxes ones. For dissipative systems, 
like the 2D Navier-Stokes equations with linear friction and stochastic forces, or the 2D Euler equations 
with linear frictions and stochastic forces, invariant measures will exist which have fluxes of conserved 
quantities. However, for two-dimensional flows, the energy flux is thought to converge to zero in the 
limit of small forcing and dissipation, by contrast to what happens for 3D flows (anomalous dissipation). 
There is thus the possibility that the limit of small forces and dissipation will be well-behaved. Some very 
interesting recent mathematical results seem to give positive insight in this direction. Our hope is 
that the invariant Young measures described in this paper could be first-order solutions in an asymptotic 
expansion of flux solutions for dissipative systems. This work is a first step in this direction, by giving a 
dynamical proof for sets of Young measures and studying their perturbed dynamics. 

Acknowledgments 

This work was supported through the ANR program STATFLOW (ANR-06-JCJC-0037-01) and through 
the ANR program STATOCEAN (ANR-09-SYSC-014). 

Appendix A: Integration in the complex plane and saddle-point approx- 
imations 

In this appendix, we evaluate the asymptotic behavior, for large N, of the following integral: 
(64) Jn (E, a) = aj dz e'^^^^^^/^v [z) with /^(z) = 11 " ^ j ' 

where {\n]n>i °^ eigenvalues of the Laplacian on the domain V. 



We first study the function /tv- It has square-root singularities at z = z„ = —iXn, with 1 < n < A^. We 
consider a complex determination of the square roots, such that fj^ has a branch cut along the half-line 
z = —ix, with X > Ai (see figure [TJ. We note that each of these singularities are integrable singularities. 

We look for an asymptotic expansion (for large A^) of /jv ■ It is a classical result [25] that 

n— >oo Ijjl 

Then 

^ 1 17^1 

^-^ \„ n->oo 47r 
n=l 

and Yl'ri=i ^/^n is a convergent series. 
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-iXi 
-1X2 



ne{z) 



Figure 1: fN{z) (|64p and integration contours used in the evaluation of integral ([M]) . 



Using 



N 



71=1 



A. 



1 ^ 



ln/jv(z) = --j;in(l 

we conclude 
(65) ln/iv(z) 

In the vicinity of z = — iAi, a direct extension of this result is 

(66) 



n=l 



\V\ ( 1 

= — (i2:)lniV + C(zz) + o — 



Z IZ 



/iv(^) ~ N'-^C{iz). 



fN{-iXi + z) 



i-iz) 



-Ci{iz), 



where C is analytic in the vicinity of —iXi. 

Function {z) is analytical except along the branch cut z = —ix, with a; > Ai. In order to compute (j64p . 
we deform the integration contour as illustrated on Figure [TJ the initial and deformed contours are shown 
in gray and red, respectively. 

It is easily checked that the contribution to Jjy of the horizontal part of the contour (see Figure [TJ are 
exponentially small, for large N. Using the change of variable z = —i{x + Xi), we obtain 



Jj^{E,a) ~ ae 



+00 



-2NaEx 



■■A/iv(-i(x + Ai)), 



-2NX,aE I 

Jo 

where A/jv is the difference between the values of /at to the left and to the right of the branch cut. Using 
([65]) we get 

JN{E,a) ~ aN »^ e ^^^^i"-^ / dx eV / AC{-i{x + Xi)). 



34 



For large A^, this last integral is clearly dominated by values of x close to zero. Then, using (|66p : 

JN{E,a) ~ ca^'^TV e / dx ev ^-^^ J x 2 

where c = —2iCi{0). Finally 

^-1/2 A^ij,, -2N\iaE 
(67) JN(E,a) ~ i^AT^ 



Appendix B: Energy-enstrophy microcanonical measure 

B-1 Energy-enstrophy microcanonical measure from a Fourier decomposition 

In this appendix, we compute the entropy for the energy-enstrophy ensemble, as well as related quantities. 
The entropy is defined by 



SK{E,r2)= lim 

TV— >-oo 



N 



with 



(68) Qk,n {E, Fa) = [ fl dujj 6{£n [uj] - E)6{r2,N [lo] - F2), 

where C does not depend on the physical parameters. It depends only on N and on the geometric factors 
{A„} and can be discarded as the entropy is always defined up to an arbitrary constant. 

We start by relaxing the enstrophy constraint: the Dirac delta in enstrophy is thus replaced with a 
Boltzmann factor in the expression of ^k,n- Then, we compute, for a > , 

(69) In {E, a) = / n S{Sn[u^] " E). 

i=i 

The relation between 1^ {E,a) and ^x,n{E,T2) shall be discussed at the end of this appendix. Let us 
use a representation as an integral in the complex plane of the remaining Dirac delta function: 



AT / N 



j=l \n=l 

^ / dfcie-^^^^^n / ^ 

J— 00 •_1 J 



27r / „ 

After computing the Gaussian integrals, we get 



— — — +00 / \ — — — 
In (E, a) = ^ / dki e 2*^^!^/^ — = jt-tJn {E, a) 
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where and /tv are defined by equation , page [ 
Using result (j67p of Appendix A, we obtain 

exp [-NG {E, a) 



(70) lNiE,a) ~ Ci(iV,{A„})C2(iV,{A„},a) 

(71) with GiE,a) =2XiaE + 



In a 



where C\ depends only on and {A„} (no dependence on the physical parameters), and C2 has no expo- 
nentially large contributions for large A^ (limjv-s>oo (liiC'2) = 0). 



From the definition (|68l) it is clear that for small AT2 and AE, 0,K,NiE,T2)AT2AE is the volume of 
the part of phase space with energy comprised between E and E + AE, and enstrophy between T2 and 
T2+AT2- We also note that Poincare inequalities impose r2,Ar[u;] > 2AiiSjv[a;] and T2[uj] > 2Xi£[ljj]. Then 
from ([68]) and (f69j) we get 

r+00 

(72) lN{E,a)= dT2exp{-NaT2)nK,N{E,T2). 

J2\iE 

It is not difficult to make for nK,N{E, V2) the same type of complex plane representation and saddle point 
approximation as the one presented for Ijv in Appendix A. However, the presentation of the computation 
would be tedious as it involves two complex auxiliary variables (similar to ki above). We thus do not 
present these computations here, but we use that a large-deviation result holds: 

(73) OK,7v(i?,r2) ~ C3(Ar,{A„})C4({A4,r2,Ar)^^Pl^^M^lM + o('iy 

Then using this last expression in (j72p and performing a saddle point approximation, we conclude that the 
thermodynamic potential (I7ip of the relaxed (canonical) ensemble is related to the entropy Sk through: 

G{E,a)= min {-Sk {E,T2) + aT2] . 

V2>2\iE 

Precisely, G{E,-) is the Legendre-Fenchel transform of Sk{E,-). It is a classical result that if G has no 
singularities, then Sk can be computed from the inverse formula SKiE,T2) = mma>Q{G{E,a) + aT2} 
(see any textbook on convex analysis or [8J). Using this inversion formula we get 

(74) SK{E,r2) = ^log(T2-2XiE) + ^^. 

We note that the entropy diverges for r2 = 2XiE, the minimal accessible enstrophy r2 for a given energy 
E. This could have been expected, as only the two microscopic states co = ±\/2Eei verify the relation 
T2 = 2XiE, as can be readily seen from the Poincare inequality. 

B-2 Computation of the entropy from the Boltzmann-Gibbs entropy 

From (I18p and (|35p we the entropy of the equilibrium state p* of the energy-enstrophy microcanonical 
measure 
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a 



■k\V\ 



V 



dr / dae-^l^l'^' 



1 fa\V\ 
- In 

2 V vr 



\V\ 

In a + In 

vr 



1 

+ 2- 



The last two terms of the rhs being generic (the entropy being defined up to a constant), we retain 



(75) 



\T>\ 



'- log a. 
2 ^ 



Appendix C: Correlation of the vorticity field 

Here, we sketch the evaluation of the order of magnitude of the two-point correlation function, for the 
vorticity field obtained from the energy-enstrophy microcanonical measure. From (j30p . we see that, for the 
measure /x^^ (|26|l . the variance = (w^ — {^n)%) of w„ is of order Following the reasoning 

that lead to equations (|30p . the joint probability distribution Pn ,n,'m{'^m'^m) (for amplitude of modes Cn 
and em) can be derived from (I27p : 

(76) Pn ,n,m {u}n,U}m) ~ Cexp 

Af— )-oo 

From this expression, for n ^ m, the correlation {{uJn^^rn))N = {^n^m— {^^n) N {^m) n) n shown 
to be of order l/N"^. Now, = X]„'^ne„(r), so {{uj{r)oj{r')))N is of order Therefore, for the 

energy-enstrophy microcanonical measure fim,K (-^ going to infinity), 

(77) ((..(r)a;(r')))=0. 

We just considered the second moment of the vorticity field. However, such a result is much more general: 
vorticity values at points r and r' are actually statistically independent for the microcanonical measure, 
as could be easily shown by extending the results (I76p and (I77p . 



iVlog Fs - 2\iE 



(An-Ai) 
Ai 



(A. 



Appendix D: Evolution of the characteristic and cumulant-generating 
functionals 

In this appendix, we derive the evolution equations for the characteristic and cumulant-generating func- 
tionals. 



37 



D-1 Characteristic functional 
For the 2D Euler equations 

In order to compute the evolution equation for the characteristic functional F[X\ = (^e^ f ^('^^'^i'^)'^'^) ^ we use 
the intermediate quantity 

^ ._ g«/ A(r)w(r)dr^ 

Using the Euler equation ([1]) and an integration by parts, we obtain 

OA 



A- j dv a;(r)v(r) • VA(r). 



We then use 



(78a) (v(r)w(r)^) = j dr' G (r, r') {u: (r') u (r) A), 

(78b) but (w(r')a;(r)A) 



5^F 



5X{y)5\{y') 



so that 



(79) ^ + M/ dr' dr VA(r) • G(r, r') 



dt JJ V / V ' '5\ir)6X{r') 

For the perturbation of an equilibrium of the 2D Euler equations 

We now apply the same tools to the case of the 2D Euler equations near a dynamical equilibrium. Consider 
O a solution to the 2D Euler equations, with Q = + eoj = Vq + ev and Vq • VOq = (see section 
14X21) . 

We want to determine the evolution equation for the characteristic functional F[X] = (^e^I just 
as we did in the case of the 2D Euler equations. Defining A := e*/^(r)^{r)dr ^nd using ([50|) , page [28l we 
get 

dA f 

— = -i / A (v • VQo + Vo • Vo; + ev • Vw) A 
= i J VX- (vi7o + Vo^^ + evo;) A . 

Since 

(80) (v(r)^) = I dr' G (r, r') {u (r') A) = -i J dr' G (r, r') 

we conclude 



(81) f + / VA(r) . [g(,.. O - "o(r)^) - V„(r) 



5X{r) 



dr' dr = 0. 
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D-2 Cumulant-generating functional 



By definition, the cumulant-generating functional is H = InF. Therefore 

dH 1 OF 5H 1 6F , 1 6^F 1 5F 5F 

and 



dt F dt ' 6X{r) F(5A(r)' 6Xir)6X{r') F (5A(r)5A(r') ^2 ^A(r) ^(rO 

Hence, using ([38|) . 

dH rr^^, ^ ^, ,J 6H SH \ . ^ 

-m^'ll ^^<--' ■ '^'C^- ■■ ' [jum^) + Sam Jam) = 

For a small perturbation of a dynamical equilibrium, using (|38|) . we get 



G{r,r'){ei{ + -^^0^7^7-7^ -Vo(r)- 



6\{r)6X{r') 6X{r) 6Xir') J ^''5A(r')y "''<5A(r 



dr' dr = 0. 



Appendix E: Uniqueness of the Vlasov statistical equilibria for repulsive 
convex potentials 

In this appendix, we prove that, for a repulsive convex potential W, the microcanonical variational problem 
for the Vlasov equation has a unique solution. This has the following consequences: for a repulsive 
potential, no phase transition exists and no ensemble inequivalence exists. The argument is extremely 
simple; it is based on the concavity of Q[p] = S[p] — j3H — J dadxdpa{a)p{a,x,p). The concavity of Q 
implies the uniqueness of the critical point of Q (the uniqueness of the equilibrium in the grand canonical 
ensemble where (3 and a are the control parameters) and thus the absence of phase transitions in the 
grand-canonical ensemble. Then a classical result of convex analysis j63) (see also a simple discussion in [8j) 
implies that there is a one-to-one relation between the constraints (i?, 7 (a)) and the Lagrange multipliers 
(/3,a(o")). Then for any energy E and distribution a {a), the microcanonical variational problem has a 
unique solution and there is no phase transition either in the microcanonical ensemble. 

As S is strictly concave and J dadxdp'j{cr)p{a,x,p) is linear, in order to prove the strict concavity of Q it 
is sufficient to prove that —j3H is concave. 

It is well-known that for systems with kinetic energy = the inverse temperature /3 is positive. 

That this is necessary can be seen directly from the normalizability of the equation for the critical state 
(j63p . Systems with possible negative temperature states are the ones with bounded phase space (see for 
instance the case of the point vortex model |57|). Then, in order to prove the concavity of —j3H it is 
sufficient to prove the convexity of H (j58p , or equivalently to prove the convexity of the potential 



V[f] ~ 1^ j '^^^Pf4'[f] ~ 2 J da^dx' m(x)m(x')M^(x — x'), 

with m{x) = f dpf{x,p). We also remark that if the system were confined by some external potential, 
because this would appear as a linear term in the functional, the convexity of H would not be affected. 
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The second-order variations of the potential read 

S^V ~ 2 J dicdx' (5m(a;)5m(x')iy(a; — x'). 

Resorting to Fourier transforms, S'm{x) = f dfe^m^e*'^^, we see that convexity of V is equivalent to 
positivity of Wk, the Fourier transform of W: 

r+oo 

Wk = / dxW{x) cos{kx); 
Jo 

recalling that W is even. 

We have Wj, = Yl'^=o ^k,n, where 

2(n+l)7r 
f * 

Wk n = dx W{x) cos{kx) 



k 

"2 



, x + 2n7r\ /tt — X + 2n7r\ /tt + x + 2n7r\ /27r — x + 2n7r 

W ; - ; - ; + W 



— / dx cos(x) 
k Jo 

I [ ' dx cos(x) [ " dy \W'{y + tt)-W' (y)] 

k Jo J x+2n2i 



k I \ k I \ k I \ k 

7T — X-\-2n7T 



Now, using that W is convex, W' is increasing and thus Wk,n > 0. Hence V k Wk > and the potential 
energy functional is convex. Therefore, H is convex, and the uniqueness of the statistical equilibria is 
proved. 

The same type of arguments could be derived in dimensions d larger than 1. We treat, for example, the 
case d = 3. Assuming the interaction potential to be isotropic, W = W{r), we have 

W{k) = — — / drW(r)e*'' "" = jTr-wr / dr rW^(r) sin(A;r). 

[ZTT) J [ZTT) k Jo 

Clearly, we can use the same reasoning as before, applying it to rW{r) rather than to W{x). We thus 
conclude that, in dimension 3, if W is a repulsive isotropic potential and if rW{r) is convex, then 
V k H^(k) > 0, the potential energy functional is convex and hence, for any value of the constraints, 
the statistical equilibria arc unique and no phase transition exists. 

Let us discuss the special case of algebraic potentials W{r) = C/r" in dimension 3. We are interested in 
systems with long-range interactions (non-integrable potentials); then we suppose a < 2. The condition 
rW (r) convex is then a > 1 (for a < 1, the Fourier transform of the potential would not be defined). 
Then the preceding argument applies to values of a with 1 < a < 2, including for instance Coulomb 
potentials. 
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